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ABSTRACT 



a Context. The diffuse gas in and nearby the Milky Way plays an important role in the evolution of the entire Galaxy. It has a complex 

■ structure characterized by neutral, weakly and highly ionized gas, dust, and molecules. 

Aims. We probe this gas through the observation of its absorption lines in the high-energy spectra of background sources. 
Methods. We use high-quality spectra of AGN Mrk 509, located at high Galactic latitudes obtained with XMM-Newton, HST and 
FUSE. We use advanced absorption models consisting of photo- and coUisional-ionization. 
Q ■ Results. We constrain the column density ratios of the different phases of the interstellar medium (ISM) and measure the abundances 

5— ( ' of C, N, O, Ne, Mg, Al, Si, S, and Fe. We detect seven discrete interstellar clouds with different velocities. One is a typical low- 

j/j , velocity cloud (LVC) and three belong to the family of the intermediate-velocity clouds (IVCs) found near the Galactic disk. These 

C^ ■ four clouds show large deviation from Solar abundances in the gas phase, mostly caused by dust depletion. The other three clouds are 

ionized high-velocity clouds (HVCs) and are located either in the Galactic environment or in the Local Group halo as suggested by 
the signatures of collisional ionization. The similar abundances and ionization structure of the HVCs suggest a common location and 
_ origin: they might belong to the remainder of an extragalactic cloud which was captured by the Galaxy. 

^ . Conclusions. We have shown that combined UV / X-ray spectroscopy is a powerful tool to investigate the ISM. In common Galactic 

O^ ■ clouds, like LVCs and IVCs, the ISM shows a complex structure consisting of at least three different temperature phases. 

Q^^ ■ Key words. ISM: abundances - ISM: dust, extinction - ISM: clouds - ISM: molecules - ISM: structure - X-rays: individuals: 

l/-^ ■ Mrk 509 - X-rays: ISM 
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CS| , 1. Introduction tic medium increases the reservoir of low metallicity gas. High- 

T— I ■ velocity clouds (HVCs) play a crucial role in this process. 

l" ■ The interstellar medium (ISM) drives the evolution of the 

.^ ; Galaxy: it is enriched with heavy elements during the course hVCs contain neutral hydrogen at velocities incompatible 

JXJ . of stellar evolution, and it also provides the source of material ^jth ^ simple model of diflFerential Galactic rotation. In practice 

5^ ; for the subsequent star formation. In the spectra of background ong uses a Local Standard of R est velocity V, sr > 90kms'' t o 

. C3 . sources the ISM gives rise to reddening and absorption lines, define HVCs (for a review see IWakker & van WoerdenlFi997h . 

The ISM show s a clear multi-phase structure (for a review, see ism absorbers with 30 < Vlsr < QOkms"' are usually de- 

[Ferriere 2001). The cold phase is a blend of dust, molecules and flngd ^s intermediate-velocity clouds (IVCs). Both IVCs and 

gas below 10 K. The warm ionized gas is weakly ionized, with hVCs might originate from our Galaxy or have an external ori- 

a temperature of ~ 10 K. The hot ionized gas is characterized gjn xhey could b e debris from G alactic fountains or infalling 

by temperatures of about 10 K. Local Group gas ( Blitz et al.l[T999 l). Of course, a Galactic foun- 

The multi-phase medium plays an important role in the evo- tain origin would imply metallicities near solar, while infall from 

lution of the Galaxy. One of the parameters that affect the stellar the Local Group or the intergalactic medium (IGM) would imply 

evolution is the metallicity of the star forming regions. Stellar metal poor gas. The presence of infalling matter is required in or- 

winds and supernovae expel part of the interstellar gas out of der to maintain star formation in the Galaxy: without a substan- 

the Galactic disk, but gravity generally stops the gas from es- tial replenishing of the gas available the star formation would 

caping, such that it falls back through th e process known as stop in a period much shorter than the Hubble Time. HVCs are 

"Galactic fountain" dShapiro & Fieldiri976l) . Gas accreted from thought to contain e nough mass to sustain star formation rate 



smaller galaxies, e.g. the Magellanic Clouds, and the intergalac- (SFR) in the Galaxy (iLehner & HowkluOl ih . 
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A multiwavelength approach provides the means to analyze 
the ISM in a complete way. For instance, combined UV and X- 
ray spectroscopy provides accurate column densities and the ve- 
locity structure of the most abundant ionic species. Ionic species 
like Si n-i v, C n-iv, N v and O v i especially, are quite common 
in HVCs (ISembach et al.ll2003h . In the last decade interstellar 
absorption lines in UV spectra of background stars h ave pro- 
vided the distances of several IVCs and HVCs (see e.g. lRichten 
I2OO6I: iLehner & Howkl 1201 ll and references therein). In sum- 
mary, all the main Hi complexes are known to be Galactic. 
IVCs are at typical distances of < 2kpc, while HVCs with 
90 < Vlsr ^ 170kms"' are found between 4-13 kpc. No 
cloud is found at Vlsr > 170km s"' towards halo stars, while 
several of them are found towards AGNs, which suggests that 
these very-high velocity clouds (VHVCs) are at larger distances. 
However all are thought to be within 40 kpc (Richter 2006), 
except the Magellanic Stream which is at about 50 kpc, much 
closer than the typical distances in the Local Group halo. A plau- 
sible scenario is that VHVCs are the next generation of HVCs in- 
falling towards the Galactic disk and slowing down during their 
interaction with it (Lehner & Howk 20 11). 

The AGN Mrk 509 has been intensively studied for both 
its intrinsic spectral features and its interesting line-of-sight. 
The Galactic latitude for Mrk 509 is -30 degrees and crosses 
the halo of our Galaxy, resulting in an important contribution 
of ionized gas. F irst evidence for ionized gas was found by 
lYork et all (Il982l) . who discovered significant absorption from 
Sin, Fen and Civ in the spectra taken with the lUE satellite. 
They also observed strong Lya absorption as well as red-shifted 
lines from Ca 11 and Na 11 in optical spectra. Blade s & MortonI 
(11983) attributed these shifted lines to corotating gas in our halo. 
This was confirmed b y the detection of neutral hydrogen with 
VrsR ~ 6 km s' (JMcGee & Newtonlll986h . ISembach et al.1 
( Il995h and lSembach et aLld 19991) discovered two sets of hnes of 

Siiii-iv and Civ with large velocities, i.e. Vlsr 230 km s"' 

and -280 km s~\ wh ich were attributed to absorption by HVCs. 
ISembach et al.l (l2003l) argued that these clouds are photo-ionized 
by the extra-Galactic background, but the presence in large quan- 
tities of hot gas such as O vi cannot be explained by the photo- 
ionization models and suggests that these clouds are probably 
interacting with the hot Galactic corona or the Local Group 
medium. To explain the ionized C, Si and O column densities 
Collins et al. (2004) concluded that the HVCs have multiple 
phases. They showed that the Siiii-iv and Civ can be produced 
by a QSO photo-ionizing background, while the O vi indicates 
colhsional ionization due to the interaction with the Galactic 
corona. We n ote that the presence of molecular H2 in the LOS 
dWa kker 2006) affects the spectral region near the O vi UV line 
and its column density estimate. 

More than fifteen years after the discovery of these HVCs, 
their structure is not yet well understood. To solve these ques- 
tions improved colhsional and photo-ionization models are 
needed to constrain the ionization processes occurring in HVCs, 
as well as a multiwavelength approach which uses all the avail- 
able archival data in order to increase the number of ions de- 
tected as well as the ionization parameter range that we can sam- 
ple, and to disentangle the molecular H2 and the O vi absorption. 

This article is one of a series of papers analyzing the deep 
and broad multiwavelength cam paign on Mrk 509. Th e overview 
of the campaign is presented in lKaastraet al.l (1201 Ibh . hereafter 
paper! Here we present the analysis of the interstellar clouds 
in the LOS towards Mrk 509 through a combined UV / X-ray 
analysis of the spectra taken with the Hubb le Space Telescope 
I Cosmic Origin Spectrograph (HST/COS) dGreenet al.1 1201 2h 



and the XMM-Newton Reflection Grating Spectrometer (RGS, 
Iden Herder et al.l l200lh . Our analysis is focused on both the 
LVC, IVCs and HVCs and the dust present along the LOS to- 
wards Mrk 509. We constrain velocity, ionization and chemical 
structure of the ISM for the different ionic species. 

The paper is organized as follows. In Sect. |2] we present the 
data. In Sect. |3]we report the relevant spectral features that we 
analyze. In Sect. |4]we describe the models we use and the re- 
sults of our analysis. The discussion and the comparison with 
previous work are given in Sect. |5] Conclusions are reported in 
Sect.lSl 



2. The data 

As part of our multiwavelength campaign we observed Mrk 509 
in the far-ultraviolet wavelength band between 1155 A and 
1760 A using the Far- Ultraviolet Channel and the medium res- 
olution gratings of the Cosmic Origins Spectrograph (COS) on- 
board the Hubble Space Telescope (HST). A detailed explana- 
tion of the observing strategy, instrum ent performance, d ata re- 
duction, and calibration can be found in lKriss et al.l (|2011), here- 
after paper VI. They also present a full-scale plot of the high- 
resolution spectrum. 

To summarize the data briefly, the observations were taken 
on 2009 December 10 and 11 simultaneously w ith the Chandr a 
LETGS observations in our campaign ( Ebr ero et al.l 1201 ll) . 
Using the COS gratings G130M and G160M, we obtained total 
exposure times of 9470 s and 16452 s, respectively. With each 
grating, we used only two diff'erent grating tilts with a single 
FP-POS for each to avoid creating gaps in the spectral regions 
of interest including the AGN outflow features. The data were 
processed with the COS calibration pipeline v2.11b at STScI. 
With central wavelengths of 1309 and 1327 for G130M, 1577 
and 1589 for G160M, and two exposures at each tilt, we ob- 
tained a sufficient variety of independent placements of the spec- 
trum on the detector to identify instrumental features in our 
high signal-to-noise spectrum. Known features such as dead 
spots were excluded from the combined data; correctable fea- 
tures such as grid-wire shadows were removed via a customized 
flat-fie ld treatment described by iKriss et aP d201ll) . iKriss et all 
(I2OIII) also describe improved wavelength calibration applied 
to the data. Finally, they deconvolved the spectra using a Lucy- 
Richardson algorithm to remove the effects of the broad wings 
of the COS line-spread function dGhavamian et al.ll2009t iKrissI 
.2011.) . This correction is essential to recover the true depths of 
narrow interstellar absorption lines; comparison to prior STIS 
spectra of Mrk 509 validates the effectiveness of the deconvo- 
lution. We use the original spectrum to identify weak spectral 
features, while we use the deconvolved spectrum for measuring 
the depth and width of identified absorption lines. 

For the X-ray portion of our study we use the stacked XMM- 
Newton RGS spectrum of Mrk 509 that consists of ten ~60ks 
individual observations. The observations and d ata reduction 
for this stacked RGS spectrum are described by iKaastra et al.l 
(l2011al) . hereafter paper II. Briefly, they used the SAS 9.0 soft- 
ware package to reduce the ten individual observations. They 
then created a fluxed spectrum for each observation and stacked 
these, using RGS 1 and 2 and both spectral orders, and taking the 
effects of the XMM-Newton multi -pointing mode into account. 

In addition to the COS and RGS spectra, we also 
use an archival spectrum obtained with the Far Ultraviolet 
Spectroscopic Explorer (FUSE) for the Ovi 1032 A absorption 
line. The FUSE data are described in paper VI. 
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For our spectra l analysis we use SPEXQ version 2.03.00 
("Kaastra et alJll99 o). SPEX is a software package optimized for 
the analysis and interpretation of high-resolution UV and X-ray 
spectra. In SPEX one can simultaneously use different compo- 
nents for continuum emission as well as absorption by photo- 
ionized or collisionally-ionized gas (see the SPEX manual). We 
scale elemental abundances to the proto-Solar values recently 
recommended by Lodders & Palme (2009). We use a linear scale 
for the abundances. Throughout this paper we use SI units which 
are also the units system used in our spectral codes, adopt Icr 
errors, and use x^ statistics unless otherwise stated. We define 
high-ionization ions to be more than four times ionized, like 
O VI- VIII, while single to four times ionized ions are defined as 
low / intermediate-ionization ions, as commonly done in X-ray 
spectroscopy. 



3. Spectral features 

3.1. The COS spectrum 

Most features in the COS spectrum are foreground absorption 
lines due to the ISM of the Galaxy. We detected absorption lines 
from seven discrete ISM components (see Table 7, paper VI). 
We also found three absorption lines due to the diffuse inter- 
galactic medium in the LOS towards the AGN. In Figs. [1] |2] 
and [3] we plot the prominent interstellar lines present in the 
COS spectrum. The lines are sorted according to their ioniza- 
tion state and ionization potential. All lines are plotted in veloc- 
ity space and in units of normalized flux. The zero-velocity of 
the scale is given by the laboratory wavelength of each transi- 
tion. For consistency and easy comparison with previous work 
we then convert velocities to the Local Standard of Rest (LSR). 
The conversion from the lab oratory to the LSR scale is given by 
vlsr =VLab+ n.l6kms-' (IB lades & Mortonlll983h . 

In the COS spectrum the ISM features are resonance lines 
and encompass up to 5 ionization states, e.g. from N i to N v. 
The strongest Unes are the N i triplet around 1200 A, the O i line 
at 1302 A, and the Feii and Alii lines at 1608.4 and 1670.8 A, 
respectively. Multiple strong transitions from both Si ii and S ii 
are shown, as well as Civ and Siiv doublets. The Cii and Sim 
lines are heavily saturated. We confirm the presence of several 
discrete absorbers, seven in the case of Civ as found in pa- 
per VI, with Vlsr of about -295, -240, -125, -65, +5, +65 and 
-h90 to -1-130 kms"'. We label these absorbers with alphabetical 
letters ranging from A to G and sort them according to their Vlsr 
(Table [U. Moreover, we adopt the standard nomenclature for the 
clouds: the absorbers that produce components A, B and C are 
HVCs as their velocity modulus is greater than 100 km s"'. The 
absorbers responsible for components D, F, and G are intermedi- 
ate velocity clouds (IVC). Component E is a low- velocity cloud 
(LVC). 



3.2. The FUSE spectrum 



FUSE observed Mrk 509 in 1999 for about 52 ks dKriss et all 
I2OOOI) and in 2000 for 62.1 ks (paper VI). We use only the last 
observation because of its higher S/N ratio. We use the FUSE 
spectrum mostly for the analysis of the absorption lines due to 
the interstellar O vi (see Fig.|5]l. 



www.sron.nl/spex 



3.3. The RGS spectrum 

The XMM-Newton RGS spectrum of Mrk 509 is complex be- 
cause most of the observed features are intrinsic to the AGN and 
some of them blend with foreground absorption lines as shown 
in paper II. Absorption and emission lines intrinsic to the AGN 
can be distinguished from the ISM fines as the former are signif- 
icantly red-shifted. Fo r a complete analysis of the AGN intrinsic 
absorption we refer toDetm ers et al.l (1201 1). hereafter paper III. 
The strongest ISM features are the O i and N i absorption lines 
at 23.5 A and 3 1 .3 A, respectively. Clear evidence of highly ion- 
ized gas is provided by the O vn-vm lines at 21.6 and 19.0 A. We 
plot the individual RGS absorption edges of neutral O, Fe, and N 
together with the most prominent high-ionization lines in Fig.|7] 



4. Spectral modeling of the ISM 

Our analysis focuses on the absorption lines of the ISM along 
the LOS towards Mrk 509. These lines are narrower than the 
features intrinsic to the AGN and always narrower than 1 A. 
The only exception is the Hi Ly o- absorption line at 1215.67 A, 
which has a FWHM of 20 A (see Fig. |4]l. In order to reproduce 
both the AGN continuum and narrow and broad emission lines, 
we use the continuum emission model of paper VI. The spectral 
ranges that are strongly affected by intrinsic absorption were ex- 
cluded as the AGN analysis is carried in other paper (Kriss et 
al. 2011). For the RGS spectral modeling we used the contin- 
uum plus line emission model (Model 2) given in paper III for 
the warm absorber First we analyze the ISM in the UV and in 
the X-rays separately due to the different resolution of the in- 
struments. The different velocity components are resolved in the 
UV spectra (see Fig.[Tl|2]and[3]l, but this it is not yet possible in 
the X-ray band. 

4.1. UV spectral fits 

We have performed the UV analysis in two steps. First, we use an 
empirical model to estimate the velocity shift and dispersion of 
each absorption line and the ionic column density. This model 
provides information on the significance and kinematics of all 
the components. It further gives hints of their ionization state 
and location within the Galaxy. In the second method we use 
physical models to provide a more realistic description of the 
several foreground absorbers. These models provide important 
information on the ionization state, abundances and the sources 
responsible for heating the ISM. A final simultaneous fit to the 
UV and X-ray spectra was performed to obtain robust results by 
covering the entire range of ionization states. 

4.1 .1 . An empirical slab model for the COS-UV spectra 

The first ISM model we use consists of seven slab components, 
which are required in order to reproduce the 7 different kine- 
matic ISM absorbers observed in the COS and FUSE spec- 
fi-a. Only C iv clearly shows all seven components (see Fig. [3]l. 
Excluding Cii, the low-ionization absorption lines require no 
more than four components. Thus low-ionization gas is not de- 
tected for the HVCs. The slab model in SPEX calculates the 
transmission of a slab of material, where all ionic column den- 
sities can be chosen independently. This has the advantage that 
the spectrum can be fitted without any prior knowledge of the 
ionization balance. The slab model jointly fits all the lines that 
are are produced by the same ion. After an acceptable spectral fit 
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Fig. 1. HST/COS absorption lines from the cold phase (ionization states I-II). The flux is normalized to the continuum emission and 
the lines are displayed in the Local Standard of Rest system. The rest frame wavelength of each resonance line is also given. The 
dotted lines represent the average velocities of the seven cloud systems jKriss et alj|2011h . which are also labeled as A to G. 



is obtained, we can compare the observed column densities with 
those predicted from photo-ionization models. Free parameters 
in the slab model are the velocity dispersion cry, the Doppler 
velocity shift v and the ionic column density A^x. 

For some UV saturated lines there is a degeneracy between 
the velocity dispersion and column density, hence the deter- 
mined column density is uncertain. For several ions we only de- 
tect saturated lines in the UV spectrum, thus we have performed 
simultaneous fits for lines with similar ionization potential, like 
H I, N I and Mg ii (see Fig.|4]i, tying their v and o-y. All lines of the 
same ion were simultaneously fitted like the four C i lines (see 
Fig. [TJ, the five Si ii and three S ii lines (see Fig. |2]i, as well as 
the Si IV, C IV and N v doublets (see Fig.[3]l. In order to simplify 
our model and to shorten the CPU time we force the compo- 
nent A, B and C (HVCs) to have the same velocity dispersion, 
and similarly for components E, F and G (IVCs, see Table [TJ. 
A preliminary fit to the strong, resolved lines of N i, Fe n, S ii, 
C IV and S iv did give consistent values for the velocity disper- 
sion (Tv within these groups. Only component D shows a rather 
smaller o-y, thus we treat it separately. In the cases where com- 
ponent D was too weak or blended, we froze o-y to the value 
obtained for Oi, the strongest line in component D. The Nv 
1238.8 A line is affected by IGM absorption (paper VI) and the 
Ovi HVC lin es are affected by the presence of molecular H2 
( IWakkeiJl2006 ). To reproduce the IGM line we added a slab of 
H I with V - 5700 km s"' , cTy =40 km s"' and very small column 



density logA^Hi - 17.2m ^ in agreement with paper VI. Two 
Gaussians were added to model the H2 absorption at 1031. 2A. 
These Gaussians have velocity shifts of -1-3 and +60kms"', 
CTy - 8kms ', and equival ent width s of 0.06 and 0.03 A, re- 
spectively in agreement with lWakked ( 12006 ). Thus, both the N v 
and O VI column densities, as given in Table [1] are corrected by 
IGM and H2 contamination. The C 11 absorption lines appear to 
be different from the other low -ionization ions like S 11 and Si 11. 
Among the low ions, only C 11 clearly shows high-velocity ab- 
sorption. Moreover, a separate fit of the C n and C i absorption 
lines provides a column density ratio C n/C i ~ 100 and a veloc- 
ity dispersion ratio crcu/o'ci > 3, which are inconsistent with 
the other column density ratios of low ions like O i-ii and N i-ii. 
This suggests that the bulk of the neutral carbon in the cold phase 
is depleted into dust grains and that most of the C 11 is provided 
by the warm phase as confirmed by our physical models (see 
Sect. 14.2.21 ) For these reasons we prefer to fit the C 11 absorption 
lines together with the intermediate ions. We empirically model 
the Hi absorption assuming that it is all in the cold phase. Our 
physical model confirms that the the neutral hydrogen present in 
the cold phase is indeed two orders of magnitude larger than that 
in the warm and hot phases. 

In Figs. |4]and|5]we show the data and best fit models for 
the most prominent lines of the cold, warm and hot ISM phases. 
The results of the best fit with the seven slab components are 
reported in Table [1] They are sorted according to the velocity 
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Fig. 2. HST/COS absorption lines from the cold phase (continued). Units are same as in Fig.[T] 



shift (columns) and ionization potential (rows). We also group 
the ions which have been fitted together: H i, N i and Mg ii; Si ii 
and Aln; Fen and Niii; Cn-iv, Siiii-iv and Nv. We generally 
detect low-ionization lines only for the LVC and IVCs (compo- 
nents E and D, F and G), thus four slab components are sufficient 
to model these low-intermediate velocity absorbers. However, 
we need seven slab components for modeling the mildly ionized 
gas providing the bulk of C ii and all two or more times ionized 
ions. We note that the cr^ on average increases with the ioniza- 
tion state. A comparison of the column densities of the three 
HVCs (component A to C) shows that the fastest component A 
is the least ionized. 

In order to simultaneously fit the absorption lines of the 
mildly ionized gas, we applied a few wavelength shifts to cer- 
tain ions to match both the IVC and HVC features. These shifts 
are likely due to residual errors in the COS and FUSE wave- 
length calibration. We applied no shift for both C iv and Si iv as 
their profiles perfectly match (see Fig. |5]l. The Si iii lines have 
been shifted by - 8 km s" ' , the C ii and Nvlinesby-10 km s" ' . 
We also shifted the FUSE O vi lines by -14.5 km s ' to obtain a 
match with the COS Civ HVC lines. Because the Ovi features 
are blended and affected by H2 absorption, in the spectral fits we 
prefer to fix all the O vi LOS velocities to those of the C iv and 
S II lines, which are well constrained. 

A physical analysis of the interstellar absorbers requires 
more realistic models. In the next subsection we test two differ- 
ent SPEX models to reproduce both the neutral and ionized gas 



phases: the hot and the xabs models, which reproduce collisional 
and photo-ionization respectively. 

4.1 .2. A collisionally-ionized model for the cold gas 

The hot model in SPEX calculates the transmission of a 
collisionally-ionized equilibrium (CIE) plasma. For a given tem- 
perature and set of abundances, the model calculates the ioniza- 
tion balance and then determines all the ionic column densities 
by scaling to the prescribed total hydrogen column density. At 
low temperatures this model mimics the n eutral interstellar gas 
(see SPEX manual and Kaastra et al. 2009). Free parameters in 
the hot model are the hydrogen column density A^h, the tempera- 
ture r, the velocity dispersion cr^ and shift v and the abundances. 
Following the ISM analysis of iPinto et alj (12010) , we have 
first modeled the cold gas with a low-temperature collisionally- 
ionized gas. The interstellar cold gas is responsible for the low 
ionization absorption lines: H i, Ci, N i, O i, Mg n, Al 11, Si 11, S 11, 
Ni II, and Fe 11. These Unes arise from the LVC and IVCs (com- 
ponents D-G) and no significant low-ionization absorption is de- 
tected at high outflow velocities (see Table[TJ. Thus, we modeled 
the cold gas with just four CIE components, one for the LVC and 
three for the IVCs. As previously assumed for the slab model, 
we take the same cr^ for the three reddest absorbers, E-F-G. We 
also adopted the same abundances for components D to G: the 
poor statistics and the line blending give degenerate fits that do 
not allow us to constrain the abundances of these four absorbers 
separately. For the strongest components E and F we performed 
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Fig. 3. HST/COS absorption lines from the warm phase (ionization states III-V). Units are same as in Fig.[T] 



a fit with decoupled abundances for S, Fe and N. The abundances 
of both components were consistent within the errors. 

At first we assumed proto-Solar abundances. The fit was not 
satisfactory for any Nh (xl > 10) due to the over-predicted ionic 
column densities: Ci ~ 3 x 10^" m"-, Sii ~ 4 x 10'^ m"^ and 
Mgii ~ 2.4 X 10^" m"^. The column densities as measured for 
both C I and Mg ii using the fit with the empirical slab models 
are much lower (see Table [T]) and the Si i lines are not detected 
in the spectrum. 

A fit with free abundances provides only a partial improve- 
ments and the main problems are the high Si ii/Si i and C ii/C i 
column density ratios. These ratios can be reproduced only 
through the combined effect of a high gas temperature and de- 
pletion onto dust. A satisfactory fit is reached when we assume 
that at least 80% of Mg and 90% of C and Si ai-e depleted from 
the gaseous phase into dust grains. 

The weakly ionized gas of the LVC and IVCs shows rather 
high temperatures of about 15 000-20 000 K as expected by the 
significant S ii column density and the absence of S i (which can- 
not be attributed to depletion into dust). This suggests that an im- 
portant fraction of hydrogen is ionized. Component E, the one 
at rest, provides 90% of the total A^hi, while component F is 
the second strongest but accounts for only ~ 5%. The measured 
(Tv are consistent with those estimated through the slab model. 
The Hi column densities measured by the two different mod- 
els are consistent with each other, but their total H i col umn is 
only about 75% of the value measured at 21 cm (Murp hy etal J 
Il996i) . which might be due to the high saturation of the UV line 



or most likely to the difference in beam-size between the 21 cm 
radio observation and the pencil-beam UV observation. Because 
the Oi 1302.2 A line is heavily saturated, it is difficult to mea- 
sure the oxygen abundance. Thus we freeze it to the proto-Solar 
value. For a complete analysis of the abundances we refer to 
Sect. 14.21 in which we fit the entire UV / X-ray dataset. We have 
also tested an alternative photo-ionization model for the cold gas 
and obtained similar results (see Sect. 14.3b . 

4.1 .3. A photo-ionized model for tine warm ionized gas 

To model the warm (mildly) ionized gas responsible for the C ii- 
IV and Si iii-iv absorption observed in all velocity components, 
we tried a collisional ionization model. However, it is not pos- 
sible to obtain a satisfactory fit with seven collisionally-ionized 
gas components for any of the HVCs (A to C) and LVC / IVCs 
(D to G). This is because collisional ionization results in a nar- 
row peak in ionization, while a wide range of ionization states is 
observed. An alternative solution might be multi-phase CIE gas, 
which we have tested without obtaining satisfactory results (see 
Sect. I4.3l l. or cooUng ga s out of ionization equilibrium (see e.g. 
iGnat & Sternberg||2007h . Here we use a photo-ionization model 
for the mildly ionized gas, which provides a good solution. 

For the photo-ionization modeling we use the xabs model 
in SPEX. The xabs model calculates the transmission of a slab 
of material, where all ionic column densities are linked through 
a photo-ionization model. The relevant parameter is the ioniza- 
tion parameter ^ = L/nr^, with L the source ionizing luminosity 
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Table 1. Spectral fits to the COS and FUSE spectra using the empirical model with 7 slab components. Dashes are non detections. 





AW 


B(«) 


Q(a) 


D(o) 


E(°) 


p(a) 


Q(o) 


O-(A) 


— 


— 





= 2.0<"' 




4.5 + 0.4 




y(b) 


— 


— 





-50 + 10 


-1 + 1 


62+1 


100 + 10 


C&^ 


— 


— 


— 


15.7 + 0.3 


17.6 + 0.1 


16.8+0.1 


16.2 + 0.1 


O-C" 


< 10*'" 


— 





2.0 + 0.4 




9 + 1 




,,(b) 


-290 + 5 


— 





-64+1 


+7 + 3 


71+5 


131 + 10 


Oi(^' 


>16.4('" 


— 





19.2 + 0.4 


21.1+0.2 


20.7+0.1 


17.1+0.1 


o-C" 


— 


— 





= 2.0<"' 




7.5 + 0.5 




yW 


— 


— 





-64 + 5 


+8+2 


70+ 1 


100 + 5 


Hi(^' 


— 


— 





— 


24.53 + 0.02 


23.2 + 0.1 


— 


NiW 


— 


— 





16.4 + 0.4 


19.5 + 0.2 


18.6 + 0.1 


17.4 + 0.1 


Mgii(^> 


— 


— 


— 


18.1+0.1 


19.8 + 0.1 


18.6 + 0.3 


18.2 + 0.1 


^(h) 


<4(rf) 


— 





= 2.0<'^> 




13+1 




yW 


-300 + 10 


— 





-65 + 5 


5+5 


65 + 5 


131 + 10 


Siii(^> 


16.7 + 0.2 


— 





17.9 + 0.3 


19.5 + 0.4 


19.4 + 0.2 


16.8 + 0.3 


Alii<^' 


16 + 1 


— 





17.3 + 0.1 


18.2 + 0.1 


17.5+0.2 


15.5 + 0.5 


o-C" 


<7.6<''> 


— 





1.1 +0.4 




7.3 + 0.3 




yW 


-300 ± 5 


— 





-72.1 +0.3 


2.8 + 0.2 


58 + 1 


89+ 1 


Fen(^) 


> 16.7 W 


— 





17.5 + 0.1 


19.3 + 0.1 


18.47 + 0.05 


17.56 + 0.03 


Niii<"> 


— 


— 


— 


16.7 + 0.2 


17.53 + 0.02 


17.30 + 0.03 


16.93 + 0.07 


^(b) 


— 


— 





3.6 + 0.6 




6.4 + 0.7 




y(b) 


— 


— 





-71.4 + 0.3 


4.9 + 0.1 


61.7+0.5 


89 + 2 


Sii<^> 


— 


— 


— 


17.7 + 0.1 


20.6 + 0.1 


18.78 + 0.02 


17.8 + 0.1 


^(b) 




12.0 + 0.2 




2.6 + 0.6 




18 + 1 




yW 


-297.6 ± 0.2 


-244 + 1 


-124 + 3 


-69+1 


6 + 1 


66 + 2 


132 + 5 


Cii*^' 


17.62 ± 0.03 


17.0 + 0.1 


< 16.4*'" 


19.2 + 0.5 


20.00 + 0.05 


19.10 + 0.05 


17.4 + 0.1 


Sim*"' 


17.53 ± 0.02 


16.59 + 0.03 


15.8+0.2 


16.0 + 0.1 


19.01 + 0.05 


18.40 + 0.06 


16.5 + 0.1 


Siiv(^> 


17.34 + 0.01 


16.35 + 0.07 


15.6 + 0.3 


15.9 + 0.2 


17.70 + 0.01 


17.43 + 0.01 


15.9 + 0.1 


Civ<'^> 


18.26 + 0.01 


17.52 + 0.02 


16.7+0.2 


16.6 + 0.2 


18.34 + 0.01 


18.14 + 0.01 


16.9 + 0.1 


Nv*'^' 


16.7 + 0.1 


16.5 + 0.2 


<16.2<'" 


16.5 + 0.2 


17.30 + 0.07 


17.2 + 0.1 


<16.7W 


^(b) 




14 + 2 




6 + 4 




40 + 6 




yW 


= -297.6*''' 


= _244W 


= -124<^> 


= -69<^' 


= 6<^> 


= 66(^' 


= 132 <^> 


vi<^> 


17.6 + 0.1 


17.74 + 0.05 


17.5+0.1 


16.2 + 0.2 


18.25 + 0.03 


18.34 + 0.08 


17.3 + 0.1 



*"' The different slab components. 

**' Both the velocity dispersion cr and the LSR v are in units of km s" 

*'^' The ionic column densities A'x are in log (m^^) units. 

*''' 2 cr lower / upper limits, respectively. 

*"' Fixed parameters. 



between 1 - lOOORyd, n the density and r the distance from 
the ionizing source. Free parameters in the xabs model are the 
hydrogen column density A^h, the ionization parameter ^, the ve- 
locity dispersion o-y, the Doppler velocity shift v, and the abun- 
dances. We have created a model consisting of seven photo- 
ionized xabs components, one for each velocity component de- 
tected. Similar to the slab and hot models, we have coupled the 
(Ty and abundances within the two LVC / IVC and HVC groups 
(see Sect. 14. l.T] and Fig.|5]i, i.e. components A-B-C and E-F-G. 

The spectral energy distribution (SED) plays a crucial role 
in determining the ionization balance in the photo-ionized layers 



(see e.g. lChakravortv et al.l2009l) . Unfortunately, the SED which 
determines the ionization balance of the interstellar gas is not 
well known. For this reason, we have tested different SEDs on 
the seven ISM absorbers (see also Fig.|6]i: 

1. local emissivity (LE) of all the galaxies and QSOs, i.e. the 
integrated emission of galaxies and QSOs as seen in the local 
Universe at z = 0; 

2. local emissivity (at z = 0) of only QSOs; 

3 . cosmic background radiation plus X-ray background as mea- 
sured by HEAOl and BeppoSAX; 

4. interstellar field SED (entirely due to starlight); 
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5. powerlaw (PL) SEDs. 

The first three SEDs represent a p urely extragalactic ioniz- 
ing s ource and have been taken from iHaardt & Madaul (1200 iL 
12012 ). The interstellar field S EP refers to th e one specified in 
the CLOUDY "Hazy" manual (lFerlandll2005l) . For the last case 
(SED 5) we have created a grid of powerlaw SEDs (Fy oc v'') 
with slopes ranging from -4.0 to -0.1 with steps of 0.05 and 
with a low-energy cut-ofF below 1 Ryd. We calculate the ion- 
ization balance through the SPEX xabsinput tool: it receives as 
input the ionizing SED and the abundances of the gas, and de- 
termines the ionization balance using Cloudy version 08.00 (see 
the SPEX manual for more details). We have used the normaliza- 
tion for the physical SEDs as given in the literature, but we note 
that the absolute normalization of the SED does not matter while 
the SED shape does (see e.g. Sect. 6.2 in the SPEX manual). In 
order to have a sufficiently broad energy band, xabsinput adds 
a minimum flux value at energies of 10 ** and lO** Ryd, which 
does not affect the results. Because we do not know the abun- 
dances a priori, we have decided to fix them to the proto-Solar 
values of lLodders & Palme (2009). We have adopted the same 
SED and abundances for the IVC group, and similarly for the 
HVC group. This choice is suggested by the fact that large devi- 
ations in ionization balance and metallicity are mostly expected 
by comparing LVC / IVCs lying in the Galactic disk with the 
halo HVCs (see below). 
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Fig. 4. Best fit to the HST/COS absorption fines of the cold gas 
with the slab model (see also Table[TJ. 




Fig. 6. SEDs which have been used in the fit of the warm photo- 
ionized gas. Arbitrary units are used in order to compare the 
shape of the SEDs. LE SED refers to QSO + galaxies. See also 
Sect. l4T3l 



Among the four physically motivated SEDs, only the local 
emissivity (LE) SED including galaxies plus QSOs provides 
a satisfactory fit with proto-Solar abundances (xy ~ 1.3, see 
Table|2]i for all seven absorbers. The other physical SEDs do not 
provide good fits even when allowing for highly non-solar abun- 
dances. The powerlaw SEDs provide results similar to the best 
fitting LE SED. For absorption from both IVCs and HVCs the 
best fit is obtained by assuming a PL SED with slope a ~ -2, 
and proto-Solar abundances. The main difference between the 
PL and LE SEDs are the derived ^ values. The ionization pa- 
rameters are systematically higher for a LE SED (see Table |2]i. 
Component A is the least ionized of the three HVCs, component 
D is the least ionized among the four IVCs. Components E and F 
have the same ionization parameter. We find that HVCs are gen- 
erally more ionized than the IVCs. The conversion from the ion- 
ization parameter ^ = L/nr^ (as usually defined in the X-rays) 
to the standard U = ny/nu (commonly used at lower energies 
in HVC works) is not trivial and depends on t he adopted SED . 
Examples of conversion factors are provided bv lNetzen(l2008l) . 

As a next step we test for deviations from Solar abundances 
for both HVCs and IVCs. However, a complete set of accurate 
abundances can only be obtained by including the X-ray spec- 
trum. Therefore we first discuss the components of the X-ray 
spectral models, before presenting a model that fits simultane- 
ously FUSE, COS and RGS data and the abundances. We note 
that the error bars in Table |2] are very small, especially for the 
column densities, due to the quality of the data. However, pos- 
sible systematic errors, like those due to deviations from proto- 
Solar abundances, are not included (see e.g. the final complete 
model in Table |5]). 



4.2. X-ray spectral fits 

The RGS spectrum of Mrk 509 is shown in detail in paper III. 
Paper II lists the strongest interstellar lines in this spectrum. In 
our analysis of this spectrum we adopt the AGN continuum and 
emission line model given in paper III, as well as the outflow 
slab model in order to subtract the absorption intrinsic to the 
AGN. The RGS X-ray spectrum is complementary to the COS- 
FUSE spectra, providing the column densities of both the weakly 
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Fig. 5. Simultaneous fit to the HST/COS and FUSE warm and hot gas components with the slab model (see also Table[T]i. 



Table 2. Best fit results for the COS spectrum using 7 photo-ionized absorbers with proto-Solar abundances. 



SED type 


A'") 


B(«) 


C *°' 


D(fl) 


E(«) 


p(o) 


QW 


VC,V<« 


= -297.6 


= -244 


= -124 


= -69 


= 6 


= 66 


= 132 




cr/") 




11.1±0.1 




10.0 ±0.5 




19.0 ±0.1 




PLtn 


log^M 


0.44 ±0.01 
-0.56 ±0.01 


0.073 ± 0.002 
-0.51 ±0.01 


0.012 ± 0.002 
-0.51 ±0.03 


0.030 ± 0.002 
-3.00 ± 0.03 


3.32 ± 0.02 
-1.70 ±0.01 


2.05 ± 0.02 
-1.70 ±0.01 


0.056 ± 0.002 
-0.56 ± 0.03 




«SED "^ 


-2.3 


-2.3 


-2.3 


-2.1 


-2.1 


-2.1 


-2.1 




o-/"' 




9.6 ±0.1 




13.3 ±0.4 




18.8 ±0.1 




LE(/' 


^H<^' 


0.54 ± 0.01 


0.064 ± 0.002 


0.009 ± 0.001 


0.033 ±0.001 


3.61 ± 0.02 


2.39 ± 0.03 


0.064 ± 0.002 




log^C) 


-0.07 ± 0.01 


0.20 ± 0.01 


0.17 ±0.08 


-2.5 ±0.1 


-1.10 ±0.01 


-1.11 ±0.01 


0.15 ±0.02 



*"' The diiferent photo-ionized xabs components. 

**' Both the velocity dispersion cr,, and vciv are in units of km s"' . The Vciv are fixed. 

*''' The hydrogen column densities A'h are in 10^^^ rcT^ units. 

*''' The ionization parameter ^ = L/nar' is in units of 10"' Wm. 

*'^' ttsED is the best fitting slope for the PL SED (see also Sect. l4.1.3] l. 

*^' PL refers to the fit assuming a power-law SED, LE refers to the Local Emissivity SED including galaxies and QSOs. 



and mildly ionized O and Ne ions, which are important to con- 
strain the warm gas. The hot gas mostly absorbs at X-ray wave- 
lengths and can be thoroughly studied only in this energy do- 
main. Moreover, some important absorption lines from neutral 
atoms are often saturated in the UV providing only lower limits 
to their column densities. In the X-rays the same ions usually 
provide non saturated lines, but with limited velocity resolution. 



4.2.1 . An empirical slab model for the RGS spectrum 

We first determine the column densities independently of ion- 
ization balance as we also did for both COS and FUSE spectra. 
This is an important check of the column densities for those ions 
with saturated lines in the UV band, in particular O i and N i. It 
also gives column densities for ions absent in the UV spectra. 
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Table 3. Empirical model for the RGS spectrum. 





NaINe'-"^ 


Nb/Ne^"> 


Nc/Ne <°' 


NdINe^"'' 


NeINe^"^ 


Ne/Ne^"^ 


NgINe^"^ 


V (kms"') 


= -298 


= -244 


= -124 


= -69 


= 6 


= 66 


= 132 


cold phase 
warm phase 
hot phase 




0.6 + 0.2 

0.22 ± 0.07 




0.10 + 0.05 

0.31 ±0.06 




0.02 ± 0.01 

0.18 + 0.06 


0.009 + 0.007 
0.017 + 0.001 
0.009 + 0.005 


= 1 


0.15+0.01 
0.58 + 0.05 
1.2 + 0.3 


0.02 ± 0.01 
0.03 ± 0.01 
0.11 ±0.04 



*"' Column density ratios for tlie cold, warm and liot phases used in the empirical model fit to the RGS spectrum. We normalize the ratios to 
component E. 



Unfortunately, the seven interstellar components which are re- 
solved in the UV spectra, are one blend in the RGS spectrum. 
This yields degeneracy when fitting the RGS spectrum and we 
propose the following solution. We create an empirical model 
with seven slab components as in Sect. 14.1.11 but we freeze the 
column density ratios to those determined from the UV spectra 
for the seven velocity components for the cold, warm and hot 
phases. These column density ratios were determined from non- 
saturated UV lines. For the cold phase (Oi-ii) we use the Fen 
and S II UV lines. For the warm phase (Om-v) we use Civ and 
Si IV lines. For the hot phase we use O vi. It is thought that most 
O VI arises from the conductive layer between the warm and the 
hot gas, but our physical models predict that in the LOS towards 
Mrk 509 at least half of the O vi belongs to the hot gas and thus 
might be a possible indicator for it (see also Sect. 14.3b . The col- 
umn density ratios adopted are displayed in Table |3] The errors 
on these average ratios are estimated from the spread in the two 
column density ratios calculated; for instance, the error on the 
Ne/Ne ratio for the cold gas is given by the difference in the re- 
spective ratios provided by the Fe ii and S ii columns. In the case 
of the hot phase we have adopted just the statistical errors on the 
O VI column densities because it is the only highly ionized ion 
(five times ionized or above) present in the UV spectrum. 

The LSR velocities of the seven slab components are fixed 
to the values measured from the UV C iv lines. Moreover, we fix 
the velocity dispersion to the averages determined with UV slab 
fits (see Table [TJ. In particular, for the low-ionization absorbers 
E, F and G we adopt a value of 9 km s \ which is the average cry 
of all the weakly ionized species. We did verify that a change of 
±5 kms"' in cry (the average scatter in cry) does not affect the 
column density estimates in the RGS fits. 

All ionic column densities that are not constrained by the 
X-ray data, but have predicted X-ray continuum absorption, are 
fixed to the values calculated from the COS slab fit. Examples 
are C i-rv and Si ii-iv. We also fix the H i column densities be- 
cause H I only produces continuum absorption in the X-ray spec- 
trum and there is no way to disentangle the different kinemat- 
ics components. We split the total A^hi - 4.44 x lO^'^m"^ 
(iMurphv et al.l 1 1 996h into four components representing D-G 
obtained above for the cold phase. The results of this model are 
listed in Table |4] and discussed in Sect. 15.21 We note that only 
the N I column densities provided by the RGS fits are larger than 
those measured in the UV spectra, while the other column den- 
sities are consistent for the two different wavelength regions. 



4.2.2. A self-consistent physical model of the ISM: 
simultaneous UV / X-ray spectral fits 

We construct a realistic ISM model consisting of a multi-phase 
structure by extending the UV ISM model. Essentially, we as- 



sume that the ISM consists of seven clouds or layers with dif- 
ferent speeds but with a similar structure the LVC and IVCs 
(components D to G) have a cold phase of collisionally-ionized 
gas and molecules, a warm photo-ionized gas phase and a hot 
collisionally-ionized gas phase. The HVCs (components A to C) 
consist only of the warm and hot gas phases as suggested by 
the absence of high-velocity cold gas in the UV spectrum (see 
Table|5]l. We simultaneously fit the RGS, COS and FUSE spec- 
tra in order to get the highest possible accuracy on line strength 
and broadening. UV saturated lines are ignored when lines from 
the same ion are not saturated in the X-ray spectrum. This is the 
case for the Oi and Ni UV lines at +Okms"' and +65 kms"', 
which are saturated while their X-ray counterparts are not. 

We take into account absorption by interstellar dust with the 
SPEX amol component. The amol model calculates the trans- 
miss ion of various mol ecules, for details see Pinto et al. (20101; 
Cos tantini et al.l (|20I2|) and the SPEX manual. The model cur- 
rently takes into account the modified edge and line structure 
around the O and Si K-edge, and the Fe K / L-edges, using mea- 
sured cross sections of various compounds taken from the X-ray 
literature. 

The velocity dispersion cTy is coupled within the component 
groups A-B-C and E-F-G as previously done (see Sect. 14.1.1b . 
The cTy of component D is still a free parameter as its lines are 
clearly narrower, especially for Oi (see Fig.[T]i. The LSR veloc- 
ities are fixed to the values estimated by the C iv slab fits (see 
TabIe[TJ. As in Sect. 14.1.21 the abundances of each gas phase are 
coupled within component groups A-B-C (HVCs) and D-E-F-G 
(IVCs). For components C, D and G it is very difficult to mea- 
sure accurate independent abundances because of their weak and 
blended profiles, while this is possible for components A, B, E, 
and F. However, we also do not want to increase the complexity 
of our model unnecessarily. Thus, we separately fitted the pro- 
files of these latter four strongest components and found that the 
abundances of components A and B were in good agreement, 
although less constrained, and the same applies for components 
E and F We also decided to keep proto-Solar abundances for all 
the HVCs because a fit with free abundances does not provide 
strong constraints besides Si/C>l. Because in terms of xl the 
physical model provides results as good as those obtained with 
the empirical slab model, we prefer to maintain our choice of 
abundances. 



4.2.3. Metal depletion 

If the cold and warm phases share the same Galactic environ- 
ment, but have different heating processes, their abundances 
might be still similar However, the ratios of the column density 
estimates that we have reported in Table[T]and|4]show strong dis- 
crepancies between neutral and ionized gas column densities. In 
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Table 4. RGS spectral fits: empirical model with 7 slab absorbers. All the units are same as in Table[T] 
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Pai- 


A 


B 


C 


D 


E 


F 


G 


V 


= -298 


= -244 


= -124 


= -69 


= 6 


= 66 


= 132 


(T 


— 


— 


— 


= 2.0 




= 9 




Hi 


— 


— 


— 


= 22.2 


= 24.6 


= 23.6 


= 22.6 


Ni 


— 


— 


— 


18.3 


20.27 + 0.06 


19.5 


18.6 


Nn 


— 


— 


— 


16.7 


18.5+0.5 


17.8 


17.0 


Oi 


— 


— 


— 


19.2 


21.4 + 0.1 


20.5 


19.6 


On 


— 


— 


— 


17.9 


19.8+0.3 


19.0 


18.2 


Nei 


— 


— 


— 


18.6 


20.7+0.1 


19.9 


19.0 


Neil 


— 


— 


— 


< 18.4 


<20.4 


<19.6 


< 18.7 


Mgi 


— 


— 


— 


< 18.7 


<20.7 


<19.9 


< 19.0 


Mgii 


— 


— 


— 


< 18.7 


<20.7 


<19.9 


< 19.0 


Fei 


— 


— 


— 


17.8 


19.7+0.3 


18.9 


18.1 


Fell 


— 


— 


— 


< 17.8 


<19.8 


<19.0 


< 18.1 


0" 




= 12 




= 2.6 




= 18 




Cv 


18.9 


18.1 


17.4 


17.3 


19.1+0.3 


18.9 


17.5 


Oiu 


19.2 


18.4 


17.6 


17.6 


19.5+0.2 


19.2 


17.8 


Oiv 


19.0 


18.2 


17.5 


17.4 


19.3+0.3 


19.0 


17.6 


Nem 


19.8 


18.9 


18.1 


18.1 


20.0 + 0.2 


19.7 


18.3 


0" 




= 14 




= 6 




= 40 




Cvi 


18.5 


18.7 


18.5 


17.1 


19.2 + 0.1 


19.3 


18.2 


Nvi 


18.3 


18.5 


18.2 


16.9 


19.0 + 0.1 


19.1 


18.0 


Nvii 


18.2 


18.3 


18.1 


16.8 


18.9 + 0.3 


19.0 


17.9 


Ovi 


18.2 


18.3 


18.1 


16.8 


18.9 + 0.3 


19.0 


17.9 


Ovii 


19.4 


19.6 


19.3 


17.9 


20.14 + 0.05 


20.2 


19.1 


Oviii 


19.0 


19.1 


18.9 


17.5 


19.7+0.1 


19.8 


18.6 


Neix 


18.9 


19.0 


18.8 


17.4 


19.6 + 0.2 


19.7 


18.6 


Nex 


18.5 


18.7 


18.4 


17.0 


19.2 + 0.5 


19.3 


18.2 


Fe XVII 


16.6 


16.7 


16.5 


15.3 


17.2 + 0.2 


17.3 


16.3 



particular the C i column density is at least two orders of mag- 
nitude smaller than that for C ii, while O i-ii and N i-ii do not 
show similar trends. It is unlikely that in the cold gas the carbon 
abundance is orders of magnitude smaller than Solar. The most 
reasonable explanation is instead that most of the neutral carbon 
is depleted in dust grains or molecules. The same applies for Fe, 
Si, Al. There are indeed no detections in the UV spectrum of 
any of their neutral transitions. Moreover, the Fe column densi- 
ties are poorly constrained by the RGS slab model (see Table |4]l 
due to the weakness of the Fe L-edge (Fig.|2ll. The column densi- 
ties estimated for oxygen and nitrogen with the diagnostic model 
are instead close to those predicted for proto-Solar abundances, 
which suggests that their depletion factors should be much lower 
than those for C and Fe. We can measure the depletion factors 
by testing for the presence of CO, H2O ice, silicates and other 
molecules through the amol model adopted. The depletion of the 
cold gas phase into dust grains argues in favor of not coupling the 
abundances of the cold and warm gas for those elements which 
are expected to be involved like C, O, Mg, Al, Si, Fe, Ni. Instead 
we couple the abundances of N, Ne and S of the cold and warm 



gas be cause these elemen ts are not expected to be highly de- 
pleted JWilms et al.ll2000h . 



4.2.4. A collisionally-ionized model for the hot ionized gas 

The main difference between the UV ISM model and the final 
ISM model is the addition of seven collisionally-ionized gas hot 
components, one for each layer They will reproduce the highly- 
ionized gas that we have previously probed with the diagnos- 
tic empirical model, see Table |4] The only free parameters for 
each component are the hydrogen column density, the temper- 
ature and the velocity dispersion. The Doppler velocities are 
coupled according to the prescriptions given above. However, 
here we couple the velocity dispersion of components D-G as 
they fully blend in all the spectra. A satisfactory fit for the hot 
gas is reached by assuming proto-Solar abundances for all the 
seven hot components. Indeed, a fit with free abundances does 
not provide any significant deviation from the abundances of 
Lodders & Palme (2009). For this reason we report the results 
obtained by assuming proto-Solar abundances for the hot gas in 
Table|5] We find that most of the hot gas is at rest and originates 
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from the the slow components E and F. We obtain only upper 
limits to the column densities of the highly-ionized gas for the 
other components. We discuss these results in Sect. 15 .21 

4.3. Alternative models 

In this section we report several tests to check the validity and 
uniqueness of our model. 

The first issue concerns the physical state of the cold weakly 
ionized gas. We have used a collisionally-ionized model in 
SPEX because at low temperatures this mimics well the cold 
interstellar gas (see Sect. l4.1.2l i. However, at these temperatures 
photo-ionization could also provide a significant contribution. 
Thus, we tried an alternative model in which four additional 
photo-ionized xabs components (E to G) substitute for the cold 
collisionally-ionized ones. The fit slightly worsens and absolute 
abundances, i.e. relative to hydrogen, are not well constrained. 
They show systematic deviations of about 25-50% from the val- 
ues obtained with the collisionally-ionized model, while abun- 
dances relative to oxygen or other well constrained ions have 
smaller deviations. In summary, a photo-ionized plasma pro- 
vides a worse fit to the cold phase, but in principle part of the 
cold component might be photo-ionized. 

Excluding C ii, ions like O i-ii and O iii-iv have been put in 
two different phases in the slab model (see Sect. 14.1. 11 1. Although 
the physical model takes into account contributions to each ion 
from any phase, we have to justify the choice adopted for the em- 
pirical slab model. At first we chose to fit the C ii together with 
the Civ and Siiii-iv because it was the only 'low' ion clearly 
showing high-velocity absorption. The physical model validates 
this choice as most of the carbon in the cold phase turns out to 
be locked into dust and the bulk of C ii comes from the warm 
photo-ionized phase. For completeness we have tested a photo- 
ionization model for each velocity component in order to fit the 
entire set of neutral to doubly ionized ions (and ignoring all 
the other lines), assuming they are from a warm weakly-ionized 
phase and decoupled from the remaining more highly ionized 
ions. However, the fit is worse than before (A^l ~ +1.0); the 
observed C ii/C i. Si ii/Si iii and S ii/S iii column density ratios are 
not well reproduced. Even if we add dust depletion for carbon 
and silicon we cannot obtain satisfactory results (at any Nh). 

We have also considered an alternative interpretation of the 
warm mildly-ionized gas that produces most of C iii-iv and Si iii- 
IV. Instead of seven photo-ionized absorbers (see Sect. 14.1.31 ) we 
used multiphase collisionally-ionized gas: each velocity com- 
ponent has been fitted with 2 or even 3 CIE components. This 
model provides the worst fit so far. It might be that the gas is out 
of equilibrium, but we cannot check this further with the current 
spectral models in SPEX. However, the goodness of fit of our 
standard model suggests that photo-ionization is a likely inter- 
pretation. 

The last important issue concerns O vi. We have so far con- 
sidered this ion, together with O vii-viii, as a tracer of the hot 
highly-ionized gas. However, as we have previously mentioned 
in Sect. 14.2.11 most of O vi is thought to arise from a cooling 
colder phase with temperatures 1 - 5 x 10^ K (see e.g. Richter 
[2006). Therefore, we have applied a collisionally-ionized model 
for this conductive phase that produces N v and O vi assuming a 
temperature of 2.1 -4.4x lO^^keV (2.5-5.1 x lO^K). However, 
this model is incompatible with another collisionally-ionized 
phase producing O vii-viii. In order to explain the O viii/O vii and 
Ne x/Ne ix ratios this highest-ionization gas also provides at least 
half of the observed O vi. In summary, either half of the O vi is 
contained in the hot ionized gas, or a significant amount of O vii 



belongs to the cooling intermediate phase. In both cases there is 
a link between the column densities of the O vi-vii ions. For this 
reason we prefer our standard model and fit together all highly- 
ionized ions including and exceeding those five times ionized, 
like O VI. 



5. Discussion 

Our analysis shows that the ISM in the LOS of Mrk 509 has a 
multi-phase structure (see Table [T]and|5]l and complex dynamics. 
The ISM structure is similar for the different velocities compo- 
nents A-G, but the LVC and IVCs show a three-phase structure 
while the HVCs show only two ionized phases. We first discuss 
the general structure of the ISM that we have probed, its con- 
stituents and chemistry. Then, we will characterize the several 
absorbers, describe their heating processes and locate them in 
the Galactic environment. 



5. 1 . ISM multi-phase structure 

In the LOS towai'ds Mrk 509 the ISM is found in different phys- 
ical and chemical forms. We found the gas at various tempera- 
tures with different heating processes. Our UV spectral model- 
ing has revealed a large sample of ionization states (see Table 
[U which is further enlarged by including the X-ray detected 
absorption. The velocity dispersion is the same within clouds 
which belong to the same family, such as the IVCs. We have in- 
deed obtained a very good fit by coupling the velocity dispersion 
of the LVC and IVC components E, F, and G (see Figs. |4] and |5]l. 
Component D is an exception. 

5.1.1. The cold phase 

We have successfully modeled the interstellar gas with three 
main phases (see Table |5]l. For the cold gas we have used a 
collisionally-ionized gas with low temperatures kT ~ 0.5 - 
2.7 eV (6 000 - 30 000 K). It provides the bulk of the neutral and 
low-ionization lines, such as O i-ii, N i-ii, Fe i-ii and S ii. We de- 
tect the cold gas only for the LVC and IVCs, no neutral gas is 
found at high speed (see Table[TJ. More than 90% of the neutral 
gas is provided by the LVC component E. Component F (Vlsr 
- 65 km s ') is the second in order of column density and has an 
ionization state consistent with the LVC (for both the cold and 
warm gas), which suggests a similar environment for the two 
clouds. 

In Sect. 14.1.21 we have already shown that the C i column 
densities are too low with respect to those of O i and N i for the 
C I / C II ratio to be explained only by heating of the cold gas 
phase. We also argued that at least 90% of C and Si from the 
cold phase is locked into molecules and dust grains. We confirm 
this with the complete physical model detailed in Table |5] We 
have indeed obtained a 2cr upper limit of 0. 1 for both the carbon 
and silicon abundances w ith respect to the proto-Solar value of 
iLodders & Palmd ( l2009h . This strongly suggests the presence of 
both carbonaceous and silicate compounds in the LOS. We have 
thus tested all the molecules available in the SPEX database and 
found the best match using a significant column of CO (note the 
feature at 23.2 A in Fig. |7]i. We get only upper limits for metal- 
lic iron, pyroxene (MgSiOj), and hematite (Fe203). We are only 
able to detect dust at rest (component E). Moreover, the broad 
profile of dust and its nearness to AGN intrinsic features makes 
the detection of weak lines due to dust at different velocities dif- 
ficult in the RGS spectrum. It is possible to estimate the CO/H2 
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Fig. 7. Mrk 509 RGS data and best fit for the complete ISM mo del (see Sect. | 4.2.2i). A ll the absorption and emission lines intrinsic to 
the AGN have been subtracted through the AGN slab model of lDetmers et alJ (l201lh . The flux is normalized to the AGN continuum 
emission. Color code of the labels: red for the hot phase, purple for the warm phase and blue for the cold phase. 



ratio in the LOS of Mrk 509 by comparing our CO column den- 
sity estimate with the H2 value measured by Wakker (2006), see 
also Table la We obtain CO/H2 ~ 0.3 (0.1 - 0.6), which is un- 
usually high for the diffuse ISM. However, we cannot be highly 
confident in this value as the CO line falls in a wavelength range 
which is strongly affected by emission and absorption features 
intrinsic to the AGN (see paper III). 



5.1.2. The warm phase 

The warm phase consists of mildly-ionized gas in photo- 
ionization equilibrium. It spreads between several ionization 
states as it provides most C ii-iv and N v. In Sect. 14.1.3! we 
tested several reasonable SEDs as photo-ionizing source for all 
the seven clouds. We obtained satisfactory fits only by taking 
into account the contribution from QSOs together with the ex- 



tragalactic emission at Z - (iHaardt & Madaul 120121) . The fit 
with the SED containing only stellar Ught is not acceptable be- 
cause this SED is not able to produce the high Civ columns, 
which we have independently measured with the slab model (see 
Table [TJ. The total local emissivity (galaxies plus quasars, LE) 
is the only SED able to photo-ionize the interstellar gas up to the 
level observed. We have also tested a diagnostic power-law (PL) 
SED and constrained ionization parameters ^ which are system- 
atically lower than those estimated with the physical LE SED 
(see Table |2]i. This might be due to the excess in softer part of 
the LE with respect to the PL SED caused by the additional ex- 
tragalactic emission. The warm gas is the interstellar component 
which is best detected at high velocities (see Fig. [3). 
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Table 5. UV / X-ray simultaneous fits with the final complete ISM model. 



Phase 


Par 


A 


B 


C 


D 


E F 


G 




v(°) 


= -297.6 


= -244 


= -124 


= -69 


= 6 =66 


= 132 




cr/") 


— 


— 


— 


<2 


9.6 ±0.1 






iVn'*' 


— 


— 


— 


0.9 ±0.1 


44 ± 1 8 ± 1 


0.09 ± 0.04 




T(eV) 


— 


— 


— 


<0.6 


1.26 ±0.01 1.24 ±0.02 


2.7 ±0.1 




C/H*'-* 


— 


— 


— 




<0.1 






N/H« 


— 


— 


— 




0.8 ±0.1 




Cold gas 


O/H^ 


— 


— 


— 




0.5 ±0.1 






Ne/H<^' 


— 


— 


— 




1.3 ±0.1 






Mg/HM 


— 


— 


— 




0.41 ± 0.04 






Al/H^ 


— 


— 


— 




0.06 ± 0.01 






Si/H<^) 


— 


— 


— 




0.14 ±0.04 






S/H<^> 


— 


— 


— 




0.67 ± 0.01 






Fe/H''-' 


— 


— 


— 




0.07 ± 0.01 






Ni/HW 


— 


— 


— 




0.07 ± 0.01 






o-,/"' 




9.8 ±0.1 




4.7 ±0.2 


18.5 ±0.1 






iVH**' 


0.52 ± 0.01 


0.065 ± 0.002 


0.010 ±0.002 


0.034 ± 0.002 


5.4 ±0.3 2.9 ±0.1 


0.013 ± 0.002 




log^<*' 


-0.03 ± 0.01 


0.17 ±0.01 


0.22 ± 0.08 


-2.3 ±0.2 


-1.13±0.01 -1.12±0.01 


0.7 ±0.2 




C/H<'-> 




= 1 






= 1 




Warm gas 


0/H''-> 
Al/H<'-> 
Si/H<'-> 
Fe/H<'-> 




= 1 






0.5 ±0.1 
0.33 ± 0.08 
0.60 ± 0.01 
0.90 ± 0.06 






cr/"' 




< 11 






16±6 




Hot gas 


^H<'" 


<0.3 


<0.4 


<0.5 


1.7 ±0.5 


19 ±3 4±1 


< 1 




T(eV) 


40 ± 10 


55 ±15 


51 ±14 


14 ±1 


160 ±10 70 ± 10 


55 ±15 




C0<'" 


— 


— 


— 


— 


7.5 ±1.5 — 


— 


Molecules 


HjOice*'" 


— 


— 


— 


— 


< 1.0 — 


— 




MgSiOa <* 


— 


— 


— 


— 


<0.8 — 


— 




MetaUicFe*'" 


— 


— 


— 


— 


<0.6 — 


— 




FezOj <'" 


— 


— 


— 


— 


<0.3 — 


— 




H2 <''•"> 


— 


— 


— 


— 


10-50 4-6 


— 



'"^ Velocity units are same as in Table [T] Doppler velocities v are fixed. 

'** The hydrogen column densities A'h are in units of 10'' m"'. The i onization parameter ^ is i n units of 10"' Wm. 

''■' Abundances ratios are in the linear proto-Solar abundance units of lLodders & Palrng l l2009h . 

*''' Molecular column densities are in units of 10^° m~ ^. 

*'"' H2 column densities are taken from lWakkej l l2006h . 



5.1.3. The hot phase 

The hot phase is characterized by highly-ionized gas with tem- 
peratures of 50 - 160 eV, i.e. 0.58 - 1.9 x lO^'K, (see Table|5l). 
It is responsible for the entire O vi-vm, N vi-vii, Ne ix and C vi 
absorption (Fig. |7|i. We detect hot gas in both HVCs, LVC and 
IVCs through the UV O vi lines (see Fig.|5j. Most of the hot gas 
is provided by the slow components E and F. The hot phases of 
the IVCs and HVCs show lower temperatures, closer to that of 
an O VI interface between the warm and the hot gas, while the 
temperature of the LVC hot phase is in agreemen t with that of 
the ty pical hot coronal gas of the Galaxy (see e.g.l Yao & WangI 



1 



l2005h . In X-rays components A-C merge with D-G and we can 
only measure upper limits for the column densities of the indi- 
vidual kinematic components (see Table |5]l. We have measured 
large O vi column densities in the FUSE spectrum (Table [1) for 
components A-C. The high column densities of ionized gas and 



the non-detection of neutral gas suggest that components A-C 
should be mostly or entkely ionized by both UV/X-ray back- 
ground and collisions with the halo of our Galaxy or the Local 
Group. 

In order to estimate the location of the hot gas of components 
E and F we perform a test on the 620 ks MOS 1-2 data (paper I). 
We selected an annular region between 10-12' around Mrk509 
in both the MOS detectors and obtained their spectra. We have 
fitted the 0.5-l.OkeV MOS spectra with a power-law continuum 
and three Gaussians to describe the O vii-viii and N vn emission 
lines. The counts and fluxes are corrected for vignetting. The line 
fluxes in photons m^^s"' sr"' are: /ovn = (4.50 + 0.35) x 10"*, 
/ovin = (l-39 + 0.15)xl04,/Nvn = (l-ll±0.12)xl0l From the 
O VIII / O VII line ratio we estimate that the hot gas has an average 
temperature of 0. 1 86 + 0.006 keV, which is close to that of com- 
ponent E (see Table|5]l. For a source with a solid angle of 1 sr, at 
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a nominal distance of 10^^ m, with the observed flux given by the 
O VIII lines, at the measured temperature above, we obtain a CIE 
emission measure Y - 1.024 x 10^" m"-*. The emission measure 
is also given by F = {ne/nu) n^^AdR, where A - 10'*'* m^ is the 
surface area, V - AdR the volume, tif^/nu = 1.198, and dR is 
the depth of the gas layer. We derive n^ dR = (8.5 + 0.9) x 10^^ 
m"^ m, which does not depend on the adopted distance. 
We consider three scenarios - In our first scenario, the gas emit- 
ting O VII and O viii is the same gas responsible for the high ion- 
ization absorption. This gas then provides the column density 
that we have measured for the absorption (Table|5]l and the emis- 
sion measure n^ dR and because A^h = "h dR we can determine 
the thickness of the gas layer dR = 1.6 x lO^pc and its den- 
sity riu = 40 m"-'. This would suggest that the hot gas is related 
to the warm-hot intergalactic medium (WHIM). So, it would be 
physically decoupled from the other two (cold and warm) ISM 
phases. In Sect. 15 .4.21 we show that components E-G can be as- 
signed to the Galactic disk. The broad and unresolved profiles of 
the high ionization O vi-viii lines are blended. As a test, we have 
fitted all the high-ionization lines with the velocities free to vary. 
In this fit only two IVCs are required (E and F) with vlsr = 30 
and 95 kms'^ kT = 65 and 1 85 eV, A^h = 9 and 6 x 10^^ m -2, re- 
spectively. Most of the O vii-viii is produced by the 185 eV gas, 
whose column density provides a depth of 700 kpc according to 
the equations mentioned above. In this case the hot gas should 
be embedded in the halo of the Local Group. This means that the 
1 .6 Mpc value determined above is not plausible. 

In our second scenario, the emitting and absorbing plasmas 
are decoupled and have different location. Generally, both the 
Local Hot Bubble (LHB), the diffuse Galactic disk/halo and ex- 
tragalactic background provide important and different contri- 
butions (see e.g. Kaastra et al. 2008; Wang 2009). The LHB is 
important at lower temperatures with a significant Ovii Is-lp 
emission line at 0.57 keV. Most of the observed O viii emission 
belongs to the hot gas of the Galactic disk and halo, which is hot- 
ter than the LHB. Extragalactic sources provide the bulk of the 
X-ray background emission above 0.7 keV. For instance, if we 
assume that the emitting gas has a density of lO'^m"^ (a com- 
mon value in the local ISM) and nemi/nabs ~ 100, we obtain that 
the emitting region has a depth of a few tens of parsec, close to 
the LHB. According to this scenario, most of the O vii emission 
originates from the Local Hot Bubble, while the bulk of the ab- 
sorption would be due to the diffuse hot interstellar gas present 
in the Galactic disk and halo. This means that our absorption 
measures are consistent with a Galactic origin for the interstellar 
clouds in the LOS towards Mrk 509 (and not W HIM) in agree- 
ment with previous UV w ork dSavage et aI.l2003i : .Sembach et alJ 
l200llComnsetani200l . 

We can probe the hot gas structure and location with a 
third alternative way. Following Yao et alj (12009 ) we predict the 
O VII-VIII line emission and column densities using their verti- 
cal exponential Galactic disk model. In their model, the gas 
density decays exponentially as a function of the height above 
the Galactic plane with a scale height of 2.8 kpc and a central 
value of 1.4x10^ m"^, and the temperature has a scale height of 
1.4kpc and a central value of 3.6x10'' K. Their model param- 
eters have large uncertainties, but can still provide useful con- 
straints. Assuming a Galactic latitude of 30 degrees, on average 
this model predicts column densities of O vii - 6^2 ^ 10'^ nT^ 
and Oviii= 2^j x lO'^m"^, originating from within a few kpc 
range from the Galactic plane. These column densities might 
contribute up to the 50% of the absorbing hot gas in our LOS 
(see TableHJi. The remaining > 50% of highly-ionized gas should 



Table 6. Column densities comparison for LVC component E. 



Par 


Method 1*''' 


Method 2<^' 


Method 3<''> 


Ni<") 


19.5 + 0.2 


20.27 ± 0.06 


20.44 + 0.06 


Oi<") 


21.1+0.2 


21.4 + 0.1 


21.05 + 0.08 


Mgii<°> 


19.8 + 0.1 


<20.7 


19.83 + 0.05 


Pen*") 


19.3 + 0.1 


< 19.8 


19.00 + 0.07 


Ovi*") 


18.25 ± 0.03 


18.9 + 0.3 


18.31+0.08 



*"' The column densities Nx are in log (m ^) units. 

<*> Method 1 is the UV slab model in Table|T] 

*'^' Method 2 refers to the X-ray slab model in Table [4] 

*''' Method 3 gives the predictions of the physical model (Table|5]l. 



belong to the more distant Galactic halo and the circumgalactic 
medium (CGM) and/or to the WHIM. 



5.2. ISM column densities 

Most of the UV and X-ray interstellar lines in our spectra are 
not trivial to disentangle. The component groups HVCs (A, B, 
C), IVCs (D, F, G) and LVC (E) provide smooth profiles, espe- 
cially in the X-ray band where their blended profiles appear like 
a single line. In UV, the HVCs are well separated from the IVCs. 
However, each component produces only a handful of strong 
and not heavily saturated lines like those of Fe n, S n and C iv 
(see Figs.[Tl|2]and|3]l. The lines in the X-ray spectrum are not 
saturated, but due to the lower resolution the different velocity 
components form one blend. Moreover, both the emission and 
absorption lines intrinsic to the AGN partly affect the interstellar 
spectral range we fitted (see paper II). 

We have computed the ionic column densities predicted by 
the three different models (see Table |6]l. We compare only the 
results obtained for component E, as this component provides 
the best constrained column densities and actually they are the 
only free parameters in the RGS fits. The reason for using an 
empirical model for both the UV and X-ray spectra is to iden- 
tify discrepancies due to saturation between the UV and X-ray 
determined ionic column densities. The UV (Method 1) and 
X-ray (Method 2) column densities are mostly consistent (see 
Table |6]). Higher values for Ni and Ovi column densities are 
obtained from the RGS spectrum. The column in UV for N i is 
clearly under-estimated because of saturation as confirmed by 
the different line ratios of the components E and F in the 1200 A 
triplet (see Figs.[T]and|4|. The O vi column density discrepancy 
is smaller, about Icr. The difference is most likely due to the 
weakness of the Ovi ls-2p line at 22 A and the spectral noise 
(see Fig.|7]i. Apart from the already discussed N i, the UV/X-ray 
model (Method 3) provides column densities consistent with the 
UV slab model. 



5.3. ISM abundances 

The ISM abundances are computed by taking into account the 
contribution from gas, dust and molecules. A thorough and ac- 
curate analysis can be done only for the rest frame component, 
which is detected in all gas and molecular phases (see Table|5]l. 
We sum the gas and dust contributions to the cold phase for com- 
ponent E, calculate the total abundances and compare them to 
those of the warm phase in Table |7] 
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Table 7. Total abundances for the cold and warm phases of IVCs. 





Cold gas 


Cold phase*** 


Warm phase 


C/H"" 


<0.1 


0.7 ±0.1 


= 1 


N/H <"■'•> 


0.8 + 0.1 


0.8 ±0.1 


0.8 ±0.1 


O/H*") 


0.5+0.1 


0.7 ±0.1 


0.5±0.1(+0.4)« 


Ne/H*"-"' 


1.3 ±0.1 


1.3 ±0.1 


1.3 ±0.1 


Mg/H*"-^' 


0.41 ± 0.04 


< 1.0 <'" 


0.4-1.0 


Al/H*"' 


0.06 ± 0.01 


<1.3<* 


0.33 ± 0.08 (+0.3)<^> 


Si/H'") 


0.14 ±0.04 


< 0.62<'" 


0.60 ± 0.01 


S/H <"•'■' 


0.67 ± 0.01 


0.67 ± 0.01 


0.67 ± 0.01 


Fe/H*"> 


0.07 ± 0.01 


< 1.0"" 


0.90 ± 0.06 


Ni/H <"•'■' 


0.07 ± 0.01 


< 1.8*'" 


0.7-1.8 



*"' Abundances ratios in proto-Solar units of iLodders & Palrng J2009h . 
**' The cold phase includes cold gas and molecules from Table[5] 
*'' Cold and warm phase abundances are coupled. 
*''' Molecules and dust predicted upper limits. 
*''' Systematic errors (see also Sect. 15. 31 ). 



There are some limitations to determining the abundances, 
for instance the non-detection of aluminates allows us only 
to give the gas contribution to the aluminum abundance. 
Aluminates absorption features are weak and a fit of the O K 
edge with additional FeAl204 just provides an upper limit to the 
aluminum abundance of about 1.3 in units of Lodders & Palme 
(|2009|) . see Table Q The CO column density can be estimated, 
while for the other molecules only upper limits are obtained. 
Thus, we can give only the predicted upper limits to the abun- 
dances of Mg, Si, and Fe. Moreover, Mg cannot be detected in 
the warm phase because most Mg ii is already provided by the 
cold gas, and intermediate-ionization ions do not provide rele- 
vant absorption lines in our energy domain, thus its abundance 
has been coupled to that of the cold gas (see also Sect. 14.2.3b . 
The same applies to Ni, for we have tested a fit with additional 
NiO molecules, which have provided an upper limit to the nickel 
abundance of about 1.8. As mentioned above the abundances of 
N, Ne and S have been coupled between the two phases as those 
elements are not suspected to be depleted from the cold phase 
into dust grains. 

A comparison between the cold and the warm phases is 
worthwhile despite these limitations. First, we find consistent 
abundances for O, Si and Fe for both phases once dust is ac- 
counted for The total iron and nitrogen abundance is close to 
the proto-Solar value. DiflFerently, a elements like O, Al, Si and S 
appear to be under-abundant. The Si and S abundance measure- 
ments have a high confidence due to the several lines that have 
been used. The Al and O abundances might be under-estimated. 
The Al column density has been measured with the 1670.8 A 
absorption line which is heavily saturated. If we fit the Al ab- 
sorption line by ignoring its saturated bottom, we get a system- 
atic eiTor on the abundance of about 0.3, i.e. about 100% (see 
Table|7ll. The oxygen low-ionization lines are highly affected by 
the AGN intrinsic features. Indeed, within Icr error of the AGN 
model (see paper III) we obtain a systematic error on the ISM 
oxygen abundance of 0.4 for the warm gas. We might also miss 
some additional contribution to the ISM oxygen from molecules 
as previously discussed (see Sect. lS.l.TT i. The carbon abundance 
in the cold phase appears to be lower than in the warm phase, 
which likely means that we miss some additional molecules like 
hydrocarbons. Models predict an optical depth of about 0.1 at 
43 A, but here the S/N ratio is too low even for the LETGS spec- 
trum in order to measure column densities (see paper V). Neon 
is the only over-abundant element, which might suggest that the 



proto-Solar Ne abundance we tise is under-estimated (see e.g. 
iDrake & Tes"tall2005Hpinto et al.ll2010l) . 

Unfortunately statistics are not high enough to estimate the 
abundances of the hot ionized gas because its lines are much 
weaker than those of the cold and warm gas. We note, how- 
ever, that the high-ionization phase can be well fitted by adopting 
proto-Solar abundances. 

5.4. The characteristics of LVC, IVCs and HVCs 

The separation between IVCs and HVCs is justified if we com- 
pare the properties of the two groups (see Table |5]l. The HVCs 
are clouds highly ionized by both UV/ X-ray background pho- 
tons and collisions with the surrounding hot gas. Neutral gas 
is absent at these high speeds in this LOS. The HVCs show 
proto-Solar Si/C ratio, while the IVCs have a carbon excess. 
Moreover, the ionization parameters ^ of the HVC warm gas are 
higher than those of the IVCs (see also Fig.|9]l. If the clouds are 
effectively photo-ionized by the same QSO / Galactic SED, con- 
sistent with our best-fit results, this might suggest that the slow 
components E and F are in a region closer to the Galactic disk, 
where less UV and X-ray photons penetrate. 

5.4.1. LVC location: metallicity method 

In order to locate the interstellar absorbers we can compare our 
measurement of the iron abundance with the Galactic metallicity 
gradient (see also Pinto et al. 2010). We use the abundance for 
the warm phase of component E as it is well constrained. We 
adopt a Galactic altitude of zero for the Sun as it is less than 
30 pc distant from the Galactic plane, which is much smaller 
than the kpc scales we are interested in. If A is the abundance of 
element X in an interstellar cloud, r the difference between the 
Galactocentric radii of the cloud and that of the Sun, and h the 
cloud altitude, then the abundance change in the line of sight can 
be written as 



dA dA , 

A A - r 1 h. 

dr dh 



(1) 



The radial and vertical dependence of the abundance can be 
expressed as Ay - AqIO"''' and A/, = AqIO"'"'' (see e.g. 
iPedicelli et al.ll2009h . where Or and o-/, are the radial and ver- 
tical sl opes of the abundance g radient, and A© is the proto-Solar 
value (iLodders & P alme 2009). We estimate the slopes by cal- 
culating the average metalli city values found in the literature 
( Yamaeata & Yoshii 1994; Maciel & Costal2009t:|Pedicelli et al.1 
2009i:iChenetal.ii2011i) . We find a, ~ -0.06 and or/, ~ -0.1 1 for 
iron, which means that the vertical gradient is steeper than the 
radial one. From Eq. ([TJ we obtain 



A 

a; 



[art 



W 



+ a,,h- 10'"''')-ln 10+1. 



(2) 



Through Eq. (|2]) we calculate the Fe abundances for a grid of 
distances d starting from 0. 1 kpc and compare them with our Fe 
abundance measurement for the warm gas of component E (see 
Table |7]l. This provides d < 0.5 kpc for component E, which 
means that the bulk of the warm interstellar absorption is local 
and belongs to the Galactic disk. The cold gas exactly follows the 
kinematics of the warm gas and thus is co-located in the Galactic 
disk as well. 
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5.4.2. LVC and IVCs location: kinematics metliod 

It is possible to probe the location of the LVC and IVCs by com- 
paring their LSR velocities with the Galactic rotation. For HVCs 
this is not possible as their LSR velocities are outside the ob- 
served range for Galactic rotation. In Fig.[8]we show the rotation 
curve of the Galaxy as function of the dis tance as measured in 
the L OS towards Mrk 509 (see chapter 9 of lBinnev & Merrifieldl 
[19981) . 

Component E, with vlsr = - 10km s' (see Table [U, is 
consistent with distances of - 1 kpc or 13 - 15 kpc. However, 
we note that the density strongly decreases with the height from 
the Galactic plane. At a distance of 14 kpc along the LOS of 
Mrk 509, the height above the Galactic plane is 7 kpc, where 
there is much less neutral gas as within a few hundred pc 
from the Galactic plane. This would imply that component E 
is thus < 1 kpc far away, consistent with our earlier conclusion 
inSect. l5AT] 

Component F shows velocities of 60 - 70 km s ' and for the 
same reason, namely the decrease in density, it must be at a dis- 
tance of 4 - 6 kpc. The large distance from both the Earth and the 
Galactic plane is confirmed by the smaller neutral column den- 
sities of component F, which are lower than those of component 
E by at least an order of magnitude (see TablelTJ. The small dis- 
tance between component E and F also explains why a very good 
fit is obtained by adopting the similar abundances, temperatures 
and ionization parameters found for both clouds (Table|5]). 

The location of the remaining two IVCs is not as easy to 
determine because the uncertainties are larger. If component D 

(vlsr 65kms"') belongs to the Galactic corotating gas, then 

it must have a distance of 25 - 30 kpc (residing at the other side 
of the Galaxy) and a height from the plane larger than 12 kpc, 
which explains why its column densities are two orders of mag- 
nitude lower than those for component E. 

It is more difficult to determine the position of compo- 
nent G. Fitting the different ions has shown a spread of about 
40km s"' in the measured velocities. Fen and Sii are consis- 
tent with Vlsr = 90kms"', while most ions show velocities 
near 130km s ' (see Table[T]i. The spread is larger than the COS 
wavelength calibration uncertainties and seems to suggest a dou- 
ble nature of component G. The smallest contribution might re- 
fer to corotating gas moving at ~ 80 - 90km s"' with a distance 
of 6 - 8 kpc. The bulk of component G is represented by the 
warm phase and moving with a velocity consistent in modulus 
with those of the HVCs. Unfortunately, the signal-to-noise ra- 
tio is small for component G and the xl does not significantly 
change by dividing this component into two subcomponents, one 
at 90 and the other at 130 km s ' . If the fast component refers to 
a Galactic fountain, then it should be close to the Galactic plane 
as these fountains are expected to reach heights of at most a few 
kpc. In fact, the high temperatures measured for the cold and 
warm gas of component G might be signatures of shocks due 
to interaction with its surrounding (see Table |5]). An alternative 
explanation for component G is that it is part of a captured ex- 
tragalactic cloud, although a satisfactory fit is obtained assuming 
the same abundances as for Galactic components D, E, and F 

5.4.3. HVCs location: ionization structure 
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Fig. 8. Rotation curve of the Galaxy in the LOS towards Mrk 509 
(see Sect. 15.4.2b . The velocities of the LVC and IVCs are also 
displayed. 
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Fig. 9. Ionization parameter versus average velocity for the ISM 
warm phase of the seven cloud systems (see also Table|5]l. 



tral gas was detected along this LOS, but'Semb ach et al.l (Il999h 
found evidence for Hi 21 cm emission at these velocities at a 
distance of about 2 degrees from Mrk 509. One plausible sce- 
nario is that the three ionized HVCs are the outskirt of one large 
captured cloud, maybe matter stripped from a satellite galaxy, 
which has fallen into the gravitational well of the Milky Way 
(see also Fig.fTOli. The hot collisionally-ionized gas suggests that 
the cloud is impacting the halo of the Galaxy or the Local Group. 
The fastest cloud, component A, is less ionized than the other 2 
HVCs, possibly because it is more distant and thus has cuiTently 
a smaller interaction with the Galactic halo. This implies that the 
interaction causes the infalhng gas to slowdown. 



A different approach is required to determine the location of the 
three HVCs (components A, B and C). All three components 
have been successfully modeled by adopting the same abun- 
dances and photo-ionizing source, extra collisional ionization, 
and a high ionization state. Interestingly, no sign of fast neu- 



5.4.4. Equilibrium in the interstellar clouds forest 

Here we propose the stability curve as an alternative method to 
probe the dynamical structure of the ISM and to test whether 
some of the components are co-located. This curve shows the 
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Fig. 10. Artistic view of the IVCs (components D-G) laying near 
the Galactic disk and HVCs (A-C) impacting the Galactic halo. 



relation between the temperature T of the gas and its pressure H, 
which is commonly defined as 



Anr^cp 6nckT 



(3) 



This curve was calculated using the SPEX xabsinput tool in 
the computation of the ionization balance of the photo-ionized 
warm phase (see Sect. 14.1.31 LE SED) and it is displayed in 
Fig, nn together with the values calculated for components A- 
G. This curve divides the T - E plane in two regions: in the 
region below the curve the heating dominates the c ooling, while 
above the curve the cooling dominates the heating (iKrolik et al.l 
Il981 b. The curve segments which have positive slopes are stable 
against thermal perturbations, while those with a negative slope 
are not. Moreover components which have the same H value on 
this curve are in pressure equilibrium and are thus likely part of 
the same interstellar structure. 

Apparently, all the seven clouds are on stable branches of 
the curve (see Fig. fTTT i. Apart from component G, which is not 
well constrained, the LVC (E), IVCs (D, F and G) and HVCs (A, 
B and C) show a rather different H value. As expected the IVC 
and HVC groups have a different nature and are not co-located. 
Components E and F share the same H and are thus supposed to 
belong to the Galactic disk environment as we have previously 
shown. Component D differs from all the others and should not 
coexist with them. As expected, HVC components A, B and C 
have quite similar H, which means that they might belong to the 
same structure. The uncertainty on the H value, due to the large 
error in the ^ value, makes the results for component G uncertain. 

5.5. Comparison with previous results 

Although the interstellar clouds show a different nature and ori- 
gin in the LOS towards Mrk 509, we have proved that the ISM 
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Fig. 11. Stability curve for the LE SED. The values for all the 
warm photo-ionized components A-G are also displayed (see 
Sects. l4T3] and l5A4l ). 



structure is well described by three phases with different ion- 
ization states (see Sect. I4.2.2I I. The cold phase is a blend of 
molecules, dust and low-ionization gas, and is not detected in 
the HVCs. The warm phase is characterized by mildly ionized 
gas in photo-ionization equilibrium and is present in both HVCs 
and IVCs. The hot phase is characterized by collisional ioniza- 
tion. The se results agree well with the current stat e of the art 
( Ferriereli200 1) as well as with what was found bv iPinto et al.l 
OOlOh in the LOS towards GS 1826-238, which is a LMXB lo- 
cated near the Galactic Center This means that the ISM of the 
Galaxy follows a certain equilibrium on Galactic scales due to 
its cooling and heating processes. 

We have measured the column densities of the different 
phases for all the seven cloud systems (see Table [TJ and com- 
pared with those found in the literature. Within the errors, our 
estimates agree we ll with the column densities found by previ- 
ous work ( Sembac h et al.ll 19951 1 19991 l20Q3t ICoUins et al.ll2()04l: 
IShull et al...2009.) . Our O vi column densities for the H VCs A-C 
and th e IVC G agree within Icr with those measured bv lFox et al] 
(2006'), who take into account the H2 absorption, which affects 
this wavelength range (see Fig. |5]l. They consider the blend of 
components A and B as one absorber, but their column density 
estimate matches the sum of the values that we have obtained 
separately for components A and B. Our total Hi column den- 
sity differs from the value measured at long wavelengths (see 
Sect. 14.1.21 ). This is most likely due to the difference between 
the beam-size of the UV and radio observations. 

Our determination of the location of the different clouds 
agrees with previous work, but improves upon the previous re- 
sults. As discussed in Sect. 15.4.31 the HVCs (A-C) are Ukely at 
the outskirts of an extragalactic cloud captured by the Galaxy 
or the Local Group a nd interacting with them as suggested by 
ISembach et al.l (Il999h . Component B has a |v| lower than com- 
ponent A maybe due to a stronger impact with the Galactic hot 
gas, as suggested by its higher ionization state. This might indi- 
cate that the external layers of the cloud slewdo wn during the in- 
fall, in agreement with the scenario proposed bv lLehner & HowM 
(1201 1ft . Such clouds, known also as VHVCs (see Sect. [1]), are 
thought to lose their H i and a significant amount of speed dur- 
ing the infall and interaction with the Galaxy. Moreover, from 
the similar total pressure and temperature properties measured 
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for the warm phase, we conclude that they likely belong to the 
same cloud structure (see Fig.fTTll. Components E (~ Okms"') 
and F (~ 60kms"') have very similar properties for their cold 
and warm phases and clearly correspond to interstellar gas which 
follows the Galactic rotation and are expected to reside at less 
than 1 and 5kpc, respectively (see Fig. [8] and Sect. | 5.3I). This 
is full y consistent with the conclusions by Blades & MortonI 
(Il983h . Component D (~ -60kms ') is mos t likely located off 
the Ga lactic plane as previously suggested bv lMorton & BladesI 
(Il986l) . but should be located at high altitudes. 



6. Conclusion 

We have presented a complete analysis of the interstellar and 
circumgalactic medium towards the AGN Mrk 509, which is 
a bright X-ray source with a high Galactic latitude of about 
-30°, through high-quality grating spectra taken with XMM- 
Newton / RGS, HST / COS and FUSE. 

- On average the ISM, as found in the form of LVCs and IVCs 
near the Galactic disk, preserves a structure consisting of 
three distinct main phases with different ionization states. 
HVCs, which usually probe the halo and CGM, have a dif- 
ferent structure and in the LOS towards Mrk 509 they show 
two main phases. 

- We have probed seven different cloud systems along the 
LOS. Our study of their kinematics and abundances has re- 
vealed the nature and location of the clouds with respect to 
the Galactic environment assuming Galactic rotation. 

- The HVCs (components A-C) are highly ionized most likely 
as the result of a larger exposure to the X-ray background and 
the interaction with the circumgalactic medium. The simi- 
lar abundances and ionization structure suggest a common 
origin for these HVCs. They might be at the outskirt of an 
extragalactic cloud captured by the Galaxy. 

- The LVC (component E) and IVC (component F) refer to in- 
terstellar gas co-rotating with the Galactic disk, as confirmed 
by the detection of dust and molecules at rest wavelengths. 

- The location of the IVCs (components D and G) is not well 
constrained, they might belong either to disk co-rotating gas 
or to Galactic fountains. 

- The column densities, ionization states and locations probed 
by our alternative models/methods agree with each other and 
with the results in the literature. All this supports our re- 
search method and justifies extension to more lines-of-sight 
and/or new observations. These analyses will provide indeed 
a better mapping of the ISM and a deeper study of its chem- 
ical composition and interaction with the entire Galaxy. 
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